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Abstract. The pair condensation energy of a finite-size superconducting particle is studied as a function of 
two control parameters. The first control parameter is the shape of the particle, and the second parameter 
is a position-dependent impurity introduced in the particle. Whereas the former parameter is known to 
induce strong fluctuations in the condensation energy, the latter control parameter is found to be a more 
gentle probe of the pairing correlations. 
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1 Introduction 

In a seminal paper, P. W. Anderson [Tj addressed the ques- 
tion of how small a metallic grain would have to be for that 
a superconducting state would cease to exist. He argued 
that quantum confinement would force the single-particle 
(sp) spectrum to become discretely resolved. The mean 
sp energy spacing will increase with decreasing size of the 
grain, until it becomes comparable to the superconduct- 
ing gap in the bulk phase. At that point, the bulk gap 
looses its significance as a clear gap between a single super- 
conducting state and a continuum of excited states, and 
the superconducting phase would evaporate. The single- 
electron transistor (SET) experiments of Ralph, Black and 
Tinkham [5] demonstrated that the energy spectrum of 
nm-scale Al particles is discretely resolved, and moreover, 
the spectrum was found to be dependent on number parity 
and externally applied magnetic fields [3l|4], establishing 
the persistence of pairing correlations at the nanoscale. 
Bardeen, Cooper and Schrieffer (BCS) [S] had identified an 
effective electron-electron pairing interaction as the driv- 
ing force behind the superconducting state in bulk mate- 
rials. A key feature of BCS theory is that the supercon- 
ducting ground state of a superconductor is modeled as 
a coherent condensation of Cooper pairs [B]. While this 
approximation is essentially valid in the bulk limit, it is 
no longer sound for finite-size systems because inaccura- 
cies induced by particle-number fluctuations become rela- 
tively large. This opened a call for theoretical approaches 
in the canonical regime, such as Lanczos diagonalisation 
|7], projected-BCS 0, Density Matrix Renormalization 
Group m, or the Richardson- Gaudin (RG) formalism [TO] . 
Richardson had shown, within the context of pairing in 
nuclear structure physics that the reduced, level- in- 



dependent, BCS Hamiltonian is exactly diagonalizable by 
means of a Bethe Ansatz product state, provided the RG 
variables, which occur as rapidities in the Ansatz, form a 
solution to the set of RG equations [T^ll3j . Later, Gaudin 
decomposed the reduced BCS Hamiltonian into a com- 
plete set of commuting conserved charges, adding it to 
the class of integrable models |14| . The main significance 
of these results is that, as long as a level- independent pair 
scattering term is considered, the BCS Hamiltonian can be 
solved exactly for a general sp spectrum, within polyno- 
mially scaling computing time. Therefore, it is a practical 
tool for the investigation of pairing correlations in meso- 
scopic systems as a function of the sp spectrum. First 
investigations were performed with a uniform sp energy 
spacing [71IH], however, studies with randomly generated 
spectra showed an enhancement of pairing correlations by 
randomness [T51[TU] . This is related to the observation that 
pairing correlations are significantly stronger around the 
Fermi level, such that a random increase of the level den- 
sity around the Fermi level will have a stronger impact 
on the mean pair correlations in a uniform sampling. The 
result of this study triggers the question whether pairing 
correlations could be enhanced in a controlled fashion. A 
sensible control parameter for the sp spectrum would be 
the shape and size of the nanoparticle. In a free- wave "par- 
ticle in a box" picture, the geometric boundary conditions 
at the surface of the nanoparticle will fix the spectrum 
of the particles. The variations in pairing correlations of 
a rectangularly shaped nanoparticle were investigated in 
this way as a function of the aspect ratio |16| . and more 
recently, the shell structure in spherically shaped nanos- 
tructures has been assessed [TTllT^rro] in connection with 
the scanning tunneling experiments (STM) on deposited 
Pb [20] and Sn [3T] nanoclusters. The reduction from three 



2 



Mario Van Raemdonck et al.: Perturbations of superconducting metallic nanoparticles 



to two dimensions, relevant for the description of pairing 
correlations in superconducting spherical coatings or mul- 
tielectron bubbles in liquid helium, has also been inves- 
tigated [22l[^ . The theoretical results in rectangular ge- 
ometries showed a strongly volatile behavior of the pairing 
condensation as a function of the shape control parameter 
(see e.g. Fig. 1 in [16]), which is understood as a direct 
consequence of rapid fluctuations in the density of active 
sp levels around the Fermi level (also referred to as the 
Debeye window) [2T| . 

In the present contribution, we will focus deeper on 
how the pairing condensation energy varies as a func- 
tion of an external control parameter. The purpose of this 
work is to investigate whether there exists such a control 
parameter which is less prone to strong fluctuations and 
allows for a more controlled manipulation of the conden- 
sation energy. Our calculations will be performed within 
the canonical RG formalism [T^[T5] . using a recently pro- 
posed pseudo-deformcd quasispin algorithm [21]. In the 
next section, we recapitulate the necessary theoretical re- 
sults of the RG formalism for metallic nanograins. The 
following section is devoted to a scrutiny of the effect of 
the fluctuating level densities on the condensation energy 
for a small and easily fathomable system. This section is 
divided into two parts. In close parallel to the work of 
Gladilin et. al. [T^], we briefly discuss geometrical effects 
on the condensation energy in a first part. In a second part, 
we introduce an impurity in otherwise clean nanograins. 

For a good review on the developments in the field of 
superconducting metallic nanograins until 2001, we would 
like to refer the reader to the review paper of von Delft 
and Ralph [25]. 

2 Richardson-Gaudin 

The reduced BCS Hamiltonian is given by 

ni m 

i7 = ^e,n.+.g ^ (1) 

i=l 'i,J = l 

with the latin indices {i, j = I . . .m} referring to a set of 
doubly-degenerate sp energies within the Debeye window 
around the Fermi level. The number operator 

fii = ajai + ata:i, (2) 

counts the number of particles within a level i and the pair 
scattering term is represented by the pair creation/anni- 
hilation operators 

Sj^ajal = (Sjy = aja., (3) 

with aj (oi) the standard fermion creation (annihilation) 
operators. The bar notation refers to the time-reversed 
partner of the corresponding operator. The set of opera- 
tors ^ and ^ span an su{2) quasispin algebra 

(4) 



with = ^{fii — 1). This algebra supports two different 
su{2) representations, corresponding to unblocked (open) 
and blocked (pair-broken) levels. In the present manuscript, 
only open levels will be considered. Richardson's result 
[T^[T^ states that the reduced BGS Hamiltonian can be 
diagonalised exactly by means of a product state of gen- 
eralised pairs, acting on the pair vacuum \6) 

N / m ct \ 

iV'(w))-n E^F^ 

a=l ' "/ 

provided the set of RG variables {x} = {xi, X2, . . . xat}, 
with N the total number of pairs, form a solution of the 
set of non-linear RG equations 

m 

l+-9E^5r^-25E^^=0' (6) 

4=1 ^ 

for a = 1, 2, . . . , A^. Since we discard blocked levels, N is 
taken as haff the total number of particles. Once the RG 
equations have been solved for {a;}, the eigenstate energy 
of the Hamiltonian is given by 

N 

It is convenient to introduce the concept of condensa- 
tion energy, which is defined as the ground state energy 
E of the system at a given pairing interaction gi, cor- 
rected by the ground state energy of the system in the 
non-interacting limit Eq 

Ec = {i^{g)\H{g)\4,{g)) - (^(0)|i7(0)|^(0)) (8) 

Because the ground state of the non-interacting system 
corresponds to a simple filling of the sp levels until the 
Fermi energy, the condensation energy of the reduced BCS 
Hamiltonian reduces to 

N N 
a—l i—1 

The benefit of using the condensation energy over the 
ground state energy is that the former quantity corrects 
for global fluctuations in the sp energy, so it is a direct 
probe for pairing correlations. 

3 Perturbations 

We will employ a simplified "particle in a box" approach 
to study the effect of perturbations. In this approach, it 
is assumed that the conductance electrons are completely 
delocalised from the atoms in the crystal, and move as 
free particles within a box, only confined by the bound- 
aries. Regarding the qualitative nature of our study, this 
approach satisfies our needs, however one should consider 
more sophisticated methods, such as Density Functional 
Theory [26j . if more realistic results are desired. 
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3.1 Geometric perturbations 

In a first part, wc study and compare the eondensation 
energy Ec within a rectangular, cylindrical and spherical 
geometry. The single-particle energies are taken as the so- 
lutions to the Schrodinger equation with the infinite-well 
sp potential V{r) 
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V{r) = 



0, Vr G 

-oo, Vr ^ 



the box 
the box 



(10) 



The wavefunctions and corresponding eigenvalues for rect- 
angular, cylindrical and spherical infinite wells can be found 
in introductory quantum mechanics textbooks |27j . For 
rectangular geometries, the sp energies are given by 



2m„ 



p 



p 



(11) 



with (nx,ny,nz) G Ng the set of quantum numbers, rrie 
the effective mass of the electron, and {lx,ly,lz) the di- 
mensions of the rectangular box. We will define a length 
scale /, such that all lengths will be given relative to /, and 
energies relative to . The sp spectrum is given in 
Figure [IJa) as a function of , with 1^ = I kept as a con- 
stant, and ly defined as such that the volume of the rectan- 
gular box also remains a constant Ixlylz ~ l^- In addition. 
Figures [TJb) and [Ijc) show the condensation energy for 
a system of A'^ = 5 pairs living in the m = 10 lowest 
sp levels of the rectangular box with respectively a weak- 
{g = -1.0[h'^/2meP]) and strong {g = -AOO^fi^ /2meP]) 
pairing interaction. Figure [IJa) not only depicts the 10 
active sp levels, but also the next 10 levels outside of the 
Debeye window (in dotted lines) to illustrate how the sp 
levels enter and leave the Debeye window as a function 
of Ix- It can be seen that the steep exits and enterings of 
the sp levels into the Debeye window give rise to strong 
fluctuations in the sp densities, especially higher up in the 
spectrum. This has an effect on the condensation energies, 
as can be inferred from Figures [TJb) andHJc). 

Before discussing the numerical results, it is worth- 
while to distill the general features in both regimes of the 
pairing strength using perturbative techniques. For the 
weak-coupling regime, standard 2nd order perturbation 
theory [27] is used to calculate the condensation energy 
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a=lb=kF+l 
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with kp the Fermi level index, and the dummy indices 
a and b running over respectively occupied and unoccu- 
pied sp levels in the non- interacting limit. For the strong- 
coupling regime, an approximate expression can be de- 
rived for the condensation energy using the RG equations 

mm 

lim Ec = gN{m-N+l)+N{{2e)-{2e)F) + 0{^), (13) 
with (2e) ~ i ^Ej the mean pair sp energy and 
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Fig. 1. The sp spectrum (a) of a rectangular box as a function 
of one of the dimensions Ix- The lowest 10 levels within the ac- 
tive Debeye window are plotted in full lines, whereas the dotted 
lines depict the next 10 levels outside the Debeye window. Fig- 
ure (b) and (c) depict the condensation energy Ec in full lines 
for a system of A'^ = 5 pairs in the m = 10 active levels of Fig- 
ure (a), for respectively the weak- and strong-couphng regime. 
The dotted lines are approximations in respective regimes (eqs. 
p2p and (|13p ). Note that the Ix axis is plotted in logarithmic 
scale to highlight the Ix -^-^ 1/lx symmetry. Deviations of this 
symmetry are due to differences in resolution on the logarith- 
mic scale. 



Fermi level. Close inspection of the functional behavior of 
expressions and ([T^ with respect to the sp spectrum 
gives away the gross features in the corresponding regime. 
In the strong-coupling regime, the condensation energy is 
dependent on (2e) and (2e) p- Therefore, the contributions 
of the sp levels on the condensation energy only depend on 
the position of the level with respect to the Fermi level, i.e. 
levels beneath the Fermi level contribute negatively with 
a weight factor (^ — 1), and levels above the Fermi level 
contribute positively with weig ht factor (^). There IS no 
direct dependency of local sp-lcvel densities on the con- 
densation energy, rather an indirect dependency entering 
via the mean sp energy above and below the Fermi level. 
For instance, it can be seen from Figure [IJc) that the 
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local increase in sp-level density around the Fermi level 
in the vicinity of Ix 0.8, 1.0 and 1.25 induces an in- 
crease of {2e)F, and therefore contributes to the pairing 
condensation. Similarly, the local increase in sp-level den- 
sity at the top of the Debeye window around Ix ~ 0.6 
and 1.6 contributes positively to (2e) and therefore de- 
creases the pairing correlations. The situation is different 
for the weak-coupling regime, where local density fluctu- 
ations around the Fermi level contribute more strongly 
than those away from the Fermi level. This can be verified 
in Figure [IJb), where, in contrast to the strong-coupling 
regime, the local density increase at Ix ~ 0.6 and 1.6 does 
not influence the condensation energy, whereas the the lo- 
cal density increase at Ix ~ 0.8, 1.0 and 1.25 considerably 
enhances the pairing correlations. 

The previous discussion highlights the importance of 
local sp-level density fluctuations entering directly, or in- 
directly into the condensation energy. Therefore, it would 
be desirable to have a more manageable control parameter 
at hand for the sp-energy levels. As is clear from Ref. [TB] 
and Figure [Ija), rectangular geometries are very prone 
to local density fluctuations and will consequently remain 
hard to control. For this reason, we performed a similar 
study of the condensation energy within a cylindrical and 
spherical configuration. Without going into much detail, 
the conclusions of these studies agreed well with the re- 
sults from the previous discussion. The sp spectrum of free 
particle waves, bounded within a cylinder with radius po 
and height 1^ is given by 



2mp 



Pi 



+ 



(14) 



the Upth root of 



with (rip, 710, riz) G N x Z x N, and a^^, 
the Bessel function J|„^|(a). It is clear that the spectrum 
(jl4p and its density fluctuations has a qualitatively sim- 
ilar dependency on the control parameter Iz as eq. (|lip . 
when volume conservation is imposed Trp^lz ~ l^- There- 
fore similar conclusions could be drawn for cylindrical as 
for rectangular geometries. 

Unfortunately, the sphere has no shape control param- 
eter if volume conservation is applied 47rpQ = l^. However, 
one can notice a gain in condensation energy in compar- 
ison with a cube and cylinder with the same dimensions 
(see Table [1]). This is nicely understood via symmetry con- 
siderations; the sphere is more symmetric than the cube, 
and will therefore exhibit more degeneracies in the sp spec- 
trum, leading to an enhancement of the pair correlations. 
Analogously, a cylinder is more symmetric than a cube 
(2D rotational vs dihedral symmetry), which generally 
translates into an enhancement for the pairing correla- 
tions, as illustrated in Table [TJ The symmetry argument 
can be tested by breaking the symmetry of the sphere 
or cylinder to spheroidal shapes. Although such a study is 
of relevance for the experiments on spherical nanodroplets 
POII^ , we do not expect significant qualitative differences 
from our study with rectangular grains and leave this sub- 
ject for future investigations. 

The conclusion of the present subsection is that the 
condensation energy of a rectangular (and cylindrical) na- 



Table 1. The condensation energy Ec (O of a cube, cylinder 
and sphere with volume l'^ for 3 different values of the pairing 
interaction strength, corresponding to a weak-, intermediate- 
and strong-couphng regime. The radius po of the cyhnder is 
fixed as such that the height — I. All calculations have been 
performed with TV = 6 pairs in the m — 12 first sp levels. 
Energies are given in units [fi^/2me/^] 



g 


cube 


cylinder 


sphere 


-1.000 


-8.851 


-6.850 


-13.067 


-10.000 


-222.254 


-226.183 


-246.037 


-20.000 


-600.819 


-608.382 


-626.128 



nograin is highly sensitive to the fluctuations in the sp- 
level densities, and that these fluctuations are rather strong 
as a function of the shape control parameter. With this 
respect, it would be interesting to find a more gentle con- 
trol parameter such that the pairing correlations can be 
probed in a more controlable fashion. In the next subsec- 
tion, we introduce impurities for this particular purpose. 



3.2 Impurities 

From a "particle in a box" perspective, an impurity can be 
modeled by means of an " obstacle" in the otherwise con- 
stant potential of the box. Let this obstacle be a Dirac S{r) 
potential. For a ID system, the potential in the Schrodinger 
equation becomes 



Vix) 



vqIS{x — xq) < X < I 

oo .T < and x > I 



(15) 



with < xo < / and vq a weighted strength parameter 
of the impurity which can be either positive or negative, 
depending on whether the impurity is considered repulsive 
or attractive. The solution to the Schrodinger equation can 
be found by solving the following transcendental equation 

^^vo sm{k[l - xo]) sin(A:xo) + kl sm{kl) = 0, (16) 
for fc G M, leading to the sp spectrum 

h^kl 



2me 



(17) 



If uq < 0, there may also exist a negative energy state, 
which is the solution of Eq. p^ . with the substitution 
ifc — > K 

Vq amh{K[l — xq]) sinh(Ka;o) + KZsinh(KZ) = 0. (18) 

In the remainder of the paper, we will only deal with re- 
pulsive impurities (vq > 0). The transcendental equation 
(|16p has a few remarkable symmetries. For instance, it 
can be verified that for xq ~ - with p G No, fc = ^ 
is always a solution of ((T6| independent from vq, as long 
as q is a multiple of p. This feature explains the quasi- 
periodic structure as a function of xq in the sp spectrum. 
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Fig. 2. In panel (a), the sp spectrum H17|) of a ID particle- in-a- 
box with a 5{x — xq) impurity at xo is depicted. Panels (b)-(d) 
show the condensation energy as a function of xo, for the weak- 
(b), intermediate- (c), and strong-coupling (d) regime. For the 
weak-, and strong-coupling regime, the approximative predic- 
tions, given in respectively Eqs. (|12l) and ()13|) . are plotted in 
dashed lines. 



which is plotted in Figure [D^a) for vq = 100.0[ft^/2TOe^^]. 
It can be seen that each level has its own frequency mod- 
ulation: the first level undergoes one full quasi-period os- 
cillation, the second level a double quasi-period oscilla- 
tion, and so on. Therefore, by probing the oscillations in 
the condensation energy, it can be inferred which sp lev- 
els contribute strongly to the final structure. We can re- 
call from the previous discussion that the condensation 



energy is more sensitive to local sp-levcl density fluctua- 
tions around the Fermi level in the weak-coupling regime, 
whereas global fluctuations contribute more in the strong- 
coupling regime. This can be observed in Figures [2Ib)-(d), 
where the condensation energy is plotted for a system 
of = 5 pairs, living in the first m = 10 sp levels of 
Figure El^a), with g = -1.0, -200.0, and -100000.0 in 
units [ft^/2meZ^]. Taking into account that the mean sp- 
energy spacing is approximately 100.0{h? /2mJ^), these 
interaction strengths correspond respectively to the weak, 
intermediate, and strong-coupling regime. In the weak- 
coupling limit (Figure [2Ib)), the condensation energy dis- 
plays 5 peaks of enhanced pairing correlations, correspond- 
ing to the 5 quasi periods of the Fermi level. On the 
other side, the 10 quasi periods at the top of the De- 
beye window are visible in the condensation energy of 
the strong-coupling limit (Figure [21Jd)). The intermediate 
regime (Figure [5Jc)) displays only 5 quasi periods, but it 
can be inferred from the shape of the modulations, that 
the signature of the top level is already present. 

The reason why the impurity is a much more gen- 
tle control parameter than the shape of a nanograin can 
be related to the relative impact of the control parame- 
ter on the sp spectrum. Whereas altering the size of the 
nanograin has a large relative effect on the available space 
of the particles-in-a-box, adding a S{x) only perturbs the 
particles marginally. The question is now whether S{r) 
perturbations are not becoming too weak when going to 
higher dimensions. In order to study this, we performed 
some exploratory calculations of the condensation energy 
with one S{r — rg) on a line (ID), in a square (2D), and 
in a cube (3D). The number of levels m = 10 and pairs 
N = 5 were chosen equal for each dimension, as well as 
the strength of the impurity vq = 100.0[ft^/2me^^], and the 
pairing strength g. These preliminary calculations pointed 
out that the condensation energy was enhanced with ap- 
proximately 20% and 25% for the 2D and 3D systems 
respectively compared to the ID case, whereas the rel- 
ative fluctuations in the condensation energy decreased 
with 75% and 85%. These numbers hint at a possible sur- 
vival of impurity induced fluctuations in the condensation 
energy in higher dimensions, but further investigations are 
required. More in particular, given the strong influence of 
the geometry on the condensation energy of the particle, 
it is unclear whether the gentle impurity-induced pertur- 
bations will be observable against the large geometric fluc- 
tuations one could encounter by picking different samples 
in a realistic experimental setting. 

The 6{x) potential has zero-range character, in con- 
trast to the spatial finite-range nature of realistic impuri- 
ties in ID systems, which may become relatively large for 
nano-sized systems [55]. Therefore, we have selected the 
following potential 



V{x) 



-k\x — xq\) < X < I 

a; < and x > I 



(19) 



to investigate the effect of the spatial extent of the impu- 
rity on the condensation energy of the ID system. Besides 
being a solution of the Helmholtz equation for screened 
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Coulomb potentials of a point-like charged particle in ID, 
the potential ([T9| acts as a distribution, including the 
6{x — xo) potential (|15p and unperturbed system in the 
/c — > oo and fc — > limit respectively. Therefore, the pa- 
rameter 1/k is a control parameter for the spatial extent 
of the impurity. We have carried out the same calcula- 
tion as in Figure [2] with the same values for vq and g, 
but for different values of k ranging from very large {S{x)- 
like) to very small (unperturbed-like) values. The results 
for k = 20/1 are depicted in Figure |3l The potential (|19p 
will deviate from the S{x — a;o) potential as k decreases, 
so the typical modulations in the sp spectrum and con- 
densation energy of the S{x — xq) case are expected to 
gradually evaporate as the potential (|19p broadens. For 
the sp spectrum, it was observed that the modulations 
were more suppressed for the high-lying states, compared 
to the low-lying states. Because the normalization of the 
potentials ([T^ and ([T5|) has been chosen equal, the po- 
tential has a finite height Vmax = volk/2, in con- 
trast with the infinite height of the 6{x) potential. As a 
result, the higher-lying excitation sp levels will be less af- 
fected by the impurity than the lower-lying sp levels (see 
Figure EJa) with Kiax = 1000[SV2me/2]). This has an 
effect on the modulations of the condensation energy in 
the strong-interaction limit. Because the modulations of 
the condensation energy in the strong-interaction limit de- 
pends approximately on the relative weighting of the sp 
levels above and beneath the Fermi level £f, the finger- 
prints of the higher sp levels will gradually disappear as k 
decreases. The value k = 20/1 has been chosen for Figure|3] 
because this is the point where the higher sp level modula- 
tions start to (visually) disappear from the condensation 
energy in the strong-interaction limit (see Figure [HJd)). 
The condensation energy in the weak-interacting limit is 
only dependent on the modulations around the Fermi level 
(see eq. ((HI)). So, we will observe modulations in the con- 
densation energy as long as the levels around the Fermi 
level are affected by the impurity. Again, this is strongly 
dependent on the relative position of the Fermi level with 
respect to the height and strength of the impurity poten- 
tial. In the limit of — > 0, all impurity induced structure 
will be lost. 



Finally, we also performed a calculation of the conden- 
sation energy for = 128 pairs living in to = 256 levels of 
the ID system with a S{x — xq) impurity. From a physics 
point of view, this particular size of system corresponds 
to a realistic number of active electron pairs within the 
Debeye window of a nanograin. The Richardson- Gaudin 
proves particularly useful in this particular situation be- 
cause the size of the Hilbert space (dim ^ 5.7 x 10''^) 
is far beyond the capabilities of standard diagonalisation 
approaches. The result of the calculation is presented in 
Figure m The effect of the larger number of particles and 
levels on the condensation energy is immediately visible 
in the modulation, which has increased to 128 quasi peri- 
ods, corresponding to the number of quasi periods of the 
Fermi level. Therefore, our analysis for the smaller system 
appears to be valid in larger systems as well. 
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Fig. 3. In panel (a), the sp spectrum (|17|l of a ID particle- 
in-a-box with a exp { — k\x — xo\) impurity at 2:0 and k = 20/Z 
is depicted. Panels (b)-(d) show the condensation energy as 
a function of xo, for the weak- (b), intermediate- (c), and 
strong-coupling (d) regime. For the weak-, and strong-coupling 
regime, the approximative predictions, given in respectively 
Eqs. (|12p and (|13p . are plotted in dashed lines. 

4 Conclusions 

In conclusion, we have studied the effect of two qualita- 
tively different control parameters on the pair condensa- 
tion energy of a finite-size superconducting particle. The 
control parameters enter into the system via the single- 
particle spectrum, which is based on a straightforward 
particle-in-a-box principle. The first control parameter is 
the shape of the particle, which induces strong fluctu- 
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Fig. 4. The condensation energy of a system consisting of 
A'^ = 128 pair in m = 256 levels, with a 5{x — xo) impurity, 
as a function of the position xq. The strength of the impu- 
rity has been chosen as vo ~ 100.0[?i^/2meZ^] and the pairing 
strength is g — — 2000.0[ft^/2me?^]. For graphical reasons, the 
condensation energy is only given in the interval xq G [0, 0.1]. 



ations in the single-particle level densities, precipitating 
into the condensation energy. By means of perturbation 
theory, it was found that the condensation energy in the 
weak-coupling regime is mainly dependent on local single- 
particle level density fluctuations, whereas the strong-coup- 
ling regime is also affected by global level density fluctua- 
tions. Introducing impurities, as a second control parame- 
ter, proved to be a more gradual probe for pairing correla- 
tions. An impurity gives a unique quasi-periodic structure 
to each single particle level as a function of the position of 
the impurity, such that it becomes possible to weigh the 
contributions of the single-particle levels to the condensa- 
tion energy by investigating the frequency of oscillations. 
Regarding the schematic nature of this study, possibilities 
to bring these results towards more realistic 2D and 3D 
system is briefly touched upon. 
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